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Abstract 

The Misner and Sharp approach to the study of gravitational col- 
lapse is extended to the viscous dissipative case in, both, the streaming 
out and the diffusion approximations. The dynamical equation is then 
coupled to causal transport equations for the heat flux, the shear and 
the bulk viscosity, in the context of Israel-Stewart theory, without ex- 
cluding the thermodynamics viscous/heat coupling coefficients. The 
result is compared with previous works where these later coefficients 
were neglected and viscosity variables were not assumed to satisfy 
causal transport equations. Prospective applications of this result to 
some astrophysical scenarios are discussed. 



1 Introduction 

The gravitational collapse of massive stars represents one of the few observ- 
able phenomena where general relativity is expected to play a relevant role. 
This fact is at the origin of the great attraction that this problem exerts on 
the comunity of the relativists, since the seminal paper by Oppenheimer and 
Snyder [1]. 

Ever since that work, much was done by researchers trying to grasp es- 
sentials aspects of this phenomenon (see [2] and references therein). However 
this endeavour proved to be difficult and uncertain. Different kind of obsta- 
cles appear, depending on the approach adopted for the modelling and/or 
on the complexity of the physical description of the fluid, assumed to form 
the selfgravitating object. All these factors in turn are conditioned by the 
relevant time scales of different physical phenomena under consideration. 

Thus, during their evolution, self-gravitating objects may pass through 
phases of intense dynamical activity, with time scales of the order of mag- 
nitude of (or even smaller than) the hydrostatic time scale, and for which 
the quasi-static approximation is clearly not reliable, e.g., the collapse of very 
massive stars [3], and the quick collapse phase preceding neutron star forma- 
tion, see for example [4] and references therein. In these cases, in which we 
are mainly concerned with here, it is mandatory to take into account terms 
which describe departure from equilibrium, i.e. a full dynamic description 
has to be used. 

We shall assume that the process is dissipative. In fact, it is already an 
established fact, that gravitational collapse is a highly dissipative process (see 
[5], [6], [7] and references therein). This dissipation is required to account 
for the very large (negative) binding energy of the resulting compact object 
(of the order of — lO^^er^f). 

Indeed, it appears that the only plausible mechanism to carry away the 
bulk of the binding energy of the collapsing star, leading to a neutron star 
or black hole is neutrino emission [8]. 

In the diffusion approximation, it is assumed that the energy flux of 
radiation (and that of thermal conduction) is proportional to the gradient of 
temperature. This assumption is in general very sensible, since the mean free 
path of particles responsible for the transfer of energy in stellar interiors is 
usually very small as compared with the typical length of the object. Thus, 
for a main sequence star as the sun, the mean free path of photons at the 



centre, is of the order of 2 cm. Also, the mean free path of trapped neutrinos 
in compact cores of densities about 10^^ g. cm. ~^ becomes smaller than the 
size of the stellar core [9, 10]. 

Furthermore, the observational data collected from supernovae 1987A 
suggests that the regime of radiation transport prevailing during the emission 
process, is closer to the diffusion approximation than to the streaming out 
limit [11]. 

Dissipative effects in the diffusion approximation are further enhanced 
by very large values of thermal conductivity, which may be of the order of 
K, ^ lO^^erg s~^ cm^^ K^^ for electron conductivity (see [12] and references 
therein) or even k ^ 10^'^erg s~^ cm,~^ K~^., for neutrino conductivity in a 
pre-supernovae event [13]. 

However in many other circumstances, the mean free path of particles 
transporting energy may be large enough as to justify the free streaming 
approximation. Therefore we shall include simultaneously both limiting cases 
of radiative transport (diffusion and streaming out), allowing for describing 
a wide range situations. 

Since we are mainly concerned with time scales that might be of the order 
of (or even smaller than) relaxation times, we have to appeal to a hyperbolic 
theory of dissipation in order to treat the transport equation for dissipative 
variables. The use of a hyperbolic theory of dissipation is further justified by 
the necessity of overcoming the difficulties inherent to parabolic theories (see 
references [5], [14]-[28] and references therein). Doing so we shall be able to 
give a description of processes occuring before thermal relaxation. 

Some years ago, Misner and Sharp [29] and Misner [30] provided a full 
account of the dynamical equations governing the adiabatic [29], and the 
dissipative relativistic collapse in the streaming out approximation [30]. 

An extension of the Misner dynamical equations as to include dissipation 
in the form of a radial heat flow (besides pure radiation), was given in [5]. 
In that work the heat flux was assumed to satisfy a causal transport equa- 
tion, but viscosity was absent and thereby the thermodynamics viscous/heat 
coupling coefficients were not taken under consideration. Furthermore, for 
simplification the fluid was assumed shearfree, in despite of the fact that the 
relevance of the shear tensor in the evolution of self-gravitaing systems has 
been brought out by many authors (see [31] and references therein). 

More recently [32], shear viscosity was introduced into the Misner ap- 
proach, but again thermodynamics viscous/heat coupling coefficients were 



neglected, and furthermore the assumed transport equation for the shear 
viscosity was the corresponding to the standard Eckart theory of relativistic 
irreversible thermodynamics [33, 34]. 

The motivation to consider viscosity effects in the study of relativistic 
gravitational collapse is well founded. In fact, though they are often excluded 
in general relativistic models of stars, they are known to play a very important 
role in the structure and evolution of neutron stars. Indeed, depending on 
the dominant process, the coefficient of shear viscosity may be as large as 
1] ^ 10^° g cm^-*^ s^^ (see [35] for a review on shear viscosity in neutron stars). 

On the other hand the coefficient of bulk viscosity may be as large as 
( ^ 10^*^ g cm~^ s~^ due to Urea processes in strange quark matter [36]. 

Similar and even larger values may be attained for two color super- 
conducting quark matter phases [37, 38] and for hybrid stars [39] (see also 
[40, 41, 42] and references therein for a review on bulk viscosity in nuclear 
and quark matter). 

The purpose of this work is to present a dynamical description of the 
gravitational collapse within the framework of the Misner approach, for the 
more general dissipative fluid distribution consistent with spherical symme- 
try. This includes the presence of both shear and bulk viscosity, with a full 
causal treatment for all dissipative variables as well as the inclusion of the 
thermodynamics viscous/heat coupling coefficients. These coefficients may 
be relevant in non-uniform stellar models [26]. 

The manuscript is organized as follows: in the next Section, besides the 
field equations, the conventions, and other useful formulae we obtain the 
resulting dynamical equation. In Section 3 transport equations in the context 
of the Miiller-Israel-Stewart theory [15, 16, 17] are obtained. The coupling of 
these equations with the dynamical equation is performed in Section 4. After 
doing that we show how the effective inertial mass density of a fluid element 
reduces by a factor which depends on dissipative variables. This result was 
already known (see [43] and references therein), but for the case in which only 
the heat flux was assumed to satisfy a causal transport equation, without the 
presence of the thermodynamics viscous/heat coupling coefficients. 

In Section 5 an expression is derived which relates the Weyl tensor with 
the density inhomogeneity and thermodynamical variables. This allows us 
to bring out the role of dissipative variables in a definition of the arrow of 
time. 

Finally the results are discussed in the last section. 
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2 The basic equations 

In this section we shall deploy the relevant equations for describing a viscous 
dissipative self-gravitating fluid. This includes a full description of the matter 
distribution, the line element, both, inside and outside the fluid boundary, 
and the field equations this line element must satisfy. 

2.1 Interior spacetime 

We consider a spherically symmetric distribution of collapsing fluid, bounded 
by a spherical surface E. We assume the fluid to undergo dissipation in the 
form of heat flow, free streaming radiation and shearing and bulk viscosity. 
Choosing comoving coordinates inside S, the general interior metric can 
be written as 

dsl = -A^de + B^dr^ + {Crf{de^ + sin^ ed<P^), (1) 

where A, B and C are functions of t and r and are assumed positive. We 
number the coordinates x^ = t, x^ = r , x"^ = 6 and x^ = 0. 

The assumed matter energy momentum T~a inside S has the form 

Tap ={f^ + P + ^) VaVp + {p + U)gap + qaVp + qpVa + djp + tt,^ (2) 

where /x is the energy density, p the pressure, 11 the bulk viscosity, g° the 
heat flux, tTq/^ the shear viscosity tensor, e the radiation density, V^ the four 
velocity of the fluid, and /° a radial null four vector. These quantities satisfy 

Ti^.V = 0, 7r[^H=0. < = 0- (3) 

In the standard irreversible thermodynamics we have [26, 44] 

Tiap = -2r](Tap, n = -C0 (4) 

where r) and ( denote the coefficient of shear and bulk viscosity, respectively, 
(Tq/3 is the shear tensor and O is the expansion. However, since we are inter- 
ested in a full causal picture of dissipative variables we shall not assume (4). 
Instead, we shall use the corresponding transport equation derived from the 
Miiller-Israel-Stewart theory. 
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The shear aap is given by 

0"a/3 = V{a;l3) + ^(aVs) - -Qha/S (5) 

where the acceleration a^ and the expansion O are given by 

aa = V^.,^V^, e = \/";« (6) 

and ha/3 = da/s + VaVf3 is the projector onto the hypersurface orthogonal to 
the four velocity . 

Since we assumed the metric (1) comoving then 

\/" = A-^6^, g" = qB-^61, P = A~^6^ + B-^6'^, (7) 

where g is a function of t and r. Also it follows from (3) that 

TToa = , nl = -2tiI = -27r33 . (8) 

In a more compact form we can write 

where x" is a unit four vector along the radial direction, satisfying 

x"xa = i, x°v; = o, x" = 5"''5r, (10) 

and r2 = |7r{. 

With (7) we obtain for (5) its non null components 

au = lB\ a,, = ^ = -\iCrfa, (11) 

6 sm d 6 

where 

and the dot stands for differentiation with respect to t, which gives the scalar 
quantity 

a,^a°^ = ^a\ (13) 



For (6) with (7) we have, 



-^■«4(l-§)' <") 



where the prime stands for r differentiation. 
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2.2 The Einstein equations 

Einstein's field equations for tlie interior spacetime (1) are given by 



G 
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(15) 



Tlie non null components of (15) with (1), (2) (7), (8) and (9) become 



Goo — SttToo 



87r(/i + e)A^ 



2^ + ^^ 
, B C C 






C" ( B' _C'\C' 2(B' _ C" 
c^Vb CJC^tAb ^C, 



^-^ 



(16) 



Gn = SttTj 



11 



iir 



p + U + e+-n 



'cr 
c 






B' = 

2 fA^ 
r I A 



B\^ 

a) 



C (C^ 



2^1 

AC 



"g 



21 



+ 



;i7) 



G. 



22 



57r 



T. 



22 



+ 



<,TX 



P + n- 



fi 



{Gt) 



~A 



B 



A" C" A' fB' 
~A^~C~^\B 



B C _AfB C\ BC' 

'b'^c^^\b^c)'^ BC_ 

C'\ B'C 1 }a' B' C'\ 
'c)~^~C^^\A~^^ ~cl 



+ 



G33 = sin^ 9G22 



(19) 



Goi = SttTqi 



-87r(g + e)AB 



fC BC CA" 



;(!-§)><-) 



observe that this last equation may be written as 



(21) 



Next, the mass function m{t,r) introduced by Misner and Sharp [29] is de- 
fined by 



{Cr 
m = — - — -n.23 



D 23 
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(22) 



2.3 The exterior spacetime and junction conditions 

Outside S we assume we have the Vaidya spacetime (i.e. we assume all out- 
going radiation is massless), described by 



ds' = - 1 



2M{v) 



dv^ - 2drdv + r^{d9^ + sin^ 



(23) 



where M{y) denote the total mass, and v is the retarded time. 

The matching of the full non-adiabatic sphere (including viscosity) to the 
Vaidya spacetime was discussed in [44]. 

From the continuity of the first and second differential forms it follows 
(see [44] for details). 



m(t,r) = M{v), 



and 
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(25) 



where = means that both sides of the equation are evaluated on E (observe 
a misprint in eq.(40) in [44] and a slight difference in notation). 
Comparing (25) with (17) and (20) one obtains 
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p + U+'^n = q. (26) 

Thus the matching of (1) and (23) on S imphes (24) and (26). 
In the context of the standard irreversible thermodynamics where (4) is 
vahd, we obtain 

P + n-i|^I,, (27) 

which reduces to eq.(41) in [44] with the appropriate change in notation. 
Observe a misprint in eq. (27) in [32] (the a appearing there is the one defined 
in [44], which is — | of the one used here and in [32]). 

2.4 Dynamical equations 

The non trivial components of the Bianchi identities , (T~"^);/3 = yield 

T;7% = -^(A + e)-i(g' + e')-2(g + e)^^ + 
'2C f ^ n\ I B / ^ 2^ 

---(^ + P + n + .--)---(^+p + n + 2. + jn 

= (28) 



and 



r;-^,„.^„ + e^.2(«'^) 



A'^ ' A BC '^ ' 

4(p' + ff + .' + |n') + l^(^ + p + n + 2. + |sj) 

To study the dynamical properties of the system, let us introduce, follow- 
ing Misner and Sharp [29], the proper time derivative D^ given by 

B.4|. (30) 
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and the proper radial derivative Dr, 



where 



D -^ ^ 



R = Cr 



(31) 



(32) 



defines the proper radius of a spherical surface inside S, as measured from 
its area. 

Using (30) we can define the velocity U of the collapsing fiuid as the 
variation of the proper radius with respect to proper time, i.e. 



U = tDtC < (in the case of collapse) . 
Then (22) can be rewritten as 



E 



{Cry 

B 



l + t/2- 



2'm{t,r) 
Cr 



1/2 



(33) 



(34) 



With (31)-(32) we can express (21) as 
47r(g + e) = E 



iflnO-.)-^; 



(35) 



Observe that in the non-dissipative, shearfree case, the equation above may 
be written, with the help of (12), (14) and (33) as 



Dn{^) = 

implying U ^ R, which describes an homologous collapse [45]. 
Next, using (16-20) and (30-32) we obtain from (22) 



(36) 



Drm = -AttR^ 



p + U + e + -n)U + {q + e)E 



and 



Drtti = A-kR^ 



/i + e+ (g + e) — 



(37) 



(38) 
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Expression (37) describes the rate of variation of the total energy inside a 
surface of radius R. On the right hand side of (37), (p + H + e + ^fi)U (in 
the case of collapse U < 0) increases the energy inside R through the rate 
of work being done by the "effective" radial pressure p + 11 + |f2 and the 
radiation pressure e. In the stationary regime where we can use the standard 
thermodynamical relation iiaiB = —2r]aai3, we recover the result obtained in 
[32]. The second term (g + e)E is the matter energy leaving the spherical 
surface. 

Equation (38) shows how the total energy enclosed varies between neigh- 
boring spherical surfaces inside the fluid distribution. The first two terms on 
the right hand side of (38), /i + e, are due to the energy density of the fluid 
element plus the energy density of the null fluid describing dissipation in the 
free streaming approximation. The last term, (g + e)U / E is negative (in the 
case of collapse) and measures the outflow of heat and radiation. 

The acceleration DtU of an infalling particle inside S can be obtained 
by using (17), (22), (30) and (34), producing 



Z,,t, = -_-4.fl(p + n + . + -f!)+— ^ 



(39) 



and then, substituting A' j A from (39) into (29), we obtain 



/i + p + n + 2e+-fi)DTf/ 



/i + p + n + 2e+-fi 
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m / 2 
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e o- 



(40) 



As it can be easily seen, the main difference between (40), and eq.( 40) 
in [32] (regarding the contributions from shear viscosity) stems from the 
Tiafi terms which now are not given by (4), but have to satisfy a transport 
equation obtained within the context of the causal dissipative theory (see 
next section). 

Thus, the factor within the round bracket on the left (which equals the 
factor on the first round bracket on the right, as it should be) represents 
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the effective inertial mass (the passive gravitational mass according to the 
equivalence principle). 

The first term on the right hand side of (40) represents the gravitational 
force. In this term, the factor within the square bracket shows how dissipation 
affects the "active" gravitational mass term. 

There are two different contributions in the second square bracket. The 
first one is just the gradient of the total "effective" pressure (which includes 
the radiation pressure and the influence of shear and bulk viscosity ). The 
second contribution comes from the local anisotropy of pressure induced by 
the radiation pressure and shear viscosity. 

The last square bracket contains different contributions due to dissipative 
processes. The third term within this bracket is positive {U < 0) showing that 
the outflow of g > and e > diminish the total energy inside the collapsing 
sphere, thereby reducing the rate of collapse. The last term describes an 
effect resulting from the coupling of the dissipative flux with the shear of the 
fluid. The effects of Dxe have been discussed in [30] and we shall not analyze 
them in detail here. 

Therefore it only remains to analyze the effects of transport equations 
when coupled to (40); we will proceed to carry on that task in the next 
section. 



3 Transport equations 

As stated in the Introduction, the main purpose of this work consists in pro- 
viding a full causal description of viscous dissipative gravitational collapse. 
This implies that all dissipative variables (noy only q) have to satisfy the 
corresponding transport equation derived from causal thermodynamics. Fur- 
thermore the thermodynamics viscous/heat coupling coefficients will not be 
neglected, as they are expected to be relevant in non-uniform stellar models 
[26]. 

Accordingly, we shall use transport equations derived from the Miiller- 
Israel-Stewart second order phenomenological theory for dissipative fluids 
[15, 16, 17]. 

This theory was proposed to overcome the pathologies [20] found in the 
approaches of Eckart [33] and Landau [34] for relativistic dissipative pro- 
cesses. The important point to retain is that this theory provides transport 



13 



equations for the dissipative variables, which are of Cattaneo type [14], lead- 
ing thereby to hyperbolic equations for dissipative perturbations. 

The starting point is the general expression for the entropy four-current, 
which in the context of the Miiller-Israel-Stewart theory, reads (see [26] for 
details) 



S^' = SnV^ 
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- iPoll^ + Piquq" + P2nu.7i'' 



\/^ aoUq^' al^x^"'q^ 
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(41) 



where n is particle number density, f3A{p,n) are thermodynamic coefficients 
for different contributions to the entropy density, and aA{p,n) are thermo- 
dynamics viscous/heat coupling coefficients. 

Next, from the Gibbs equation and Bianchi identities, it follows that 
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(42) 



Finally, by the standard procedure, the constitutive transport equations 
follow from the requirement 5'";^ > 

T,U^V'^ + U = -C& + aoCq^^-^CT(^vA H, (43) 



Tih^.^V + q. 



— K 
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(44) 
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and 
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T2h'^h''p'n^i,.pV^ + 7ra/3 = -^tjaafS + '^1l<^l<l<fS;a> " V^ I IT^^" ) ^a/3 (45) 



with 



q<f3;a> = h^pK ( o(9m;^ + (li^;f^) - ^9^;^/^'"' V 



and where the relaxational times are given by 



(46) 



ro = C/3o, ri = ft:T/3i, r^ = 27^(32. (47) 

Equations (43)-(45) reduce to equations (2.21), (2.22) and (2.23) in [26] 
when thermodynamics couphng coefficients vanish. 

In our case each of the equations (43)-(45) has only one independent 
component, they read 
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We shall now proceed to couple transport equations in the form above, to 
the dynamical equation (40), in order to bring out the effects of dissipation 
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on the dynamics of the collapsing sphere. For that purpose, let us first 
substitute (49) into (40), then we obtain, after some rearrangements. 
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(61) 



where Fgrav and Fhyd are defined by 
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(54) 



Next we express by means of (48) and feed this back into (51), obtain- 



mg: 
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16 



where A is given by 

A = «„C,(. + P + n + 2. + ?n)"(|^). (56) 

Thus, once transport equations have been taken into account, then the 
inertial energy density and the "passive gravitational mass density" , appear 
diminished by the factor 1 — A + A. This result generalizes the one obtained 
in [32] , by means of a complete causal treatment of all dissipative variables 
and the inclusion of thermodynamics viscous/heat coupling coefficients. 



4 The Weyl tensor 

In this section we shall find some interesting relationships linking the Weyl 
tensor with matter variables, from which we shall extract some conclusions 
about the arrow of time. 

From the Weyl tensor we may construct the Weyl scalar C^ = C^^"'^ C ap-yS 
which can be given in terms of the Kretchman scalar IZ = R'^'^^^Rap-yS, the 
Ricci tensor R^js and the curvature scalar R by 

C^ = 7^ - 2i?°^i?„/3 + ^R^. (57) 

With the help of the formulae given in the Appendix of [32] and the field 
equations, we may write (57) as 

£ = m R^ if^-^), (58) 

o 

where S is given by 



From (58) with (37) and (38) we have 
Dr£ = AirR 



TJ ID- 



(60) 



From (60) we obtain at once for the non-dissipative, perfect fiuid case 

4:71 

Dr£ + —R^Drix = 0, (61) 
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implying that -Dr/i = produces C = (using the regular axis condition), 
and conversely the conformally flat condition implies homogeneity in the 
energy density. 

Since tidal forces tend to make the gravitating fluid more inhomogeneous 
as the evolution proceeds, a relationship like (61) led Penrose to propose a 
gravitational arrow of time in terms of the Weyl tensor [46]. 

However the fact that such a relationship is no longer valid in the pres- 
ence of local anisotropy of the pressure and/or dissipative processes, already 
discussed in [6], explains its failure in scenarios where the above-mentioned 
factors are present [47]. 

Here we see how shear viscosity and dissipative fluxes affect the link 
between the Weyl tensor and density inhomogeneity. From the above it 
is evident that density inhomogeneities may appear in a conformally flat 
spacetime, if dissipative processes occur. Examples of this kind have been 
presented in [48]. 

5 Conclusions 

We have established the set of equations governing the structure and evolu- 
tion of self-gravitating spherically symmetric dissipative viscous fluids. 

Dissipative variables have been assumed to satisfy transport equations 
derived from causal thermodynamics, and viscous/heat coupling coefficients 
have been included. 

As a result of this aproach we obtain a dynamic equation (55) which shows 
the influence of dissipative variables and viscous/heat coupling coefficients 
on the value of the "effective" inertial mass ( "passive" gravitational mass). 

In a presupernovae event, dissipative parameters (in particular k,) may be 
large enough as to poduce a significative decreasing of the gravitational force 
term, resulting in a reversal of the collapse. A numerical model exhibiting 
this kind of "bouncing" has been presented in [49]. 

Nevertheless, the role that this effect might play in the outcome of gravi- 
tational collapse of massive stars will critically depend on specific numerical 
values of those quantities. Such estimations are, however, well beyond the 
scope of this work. 

Here we just want to display the way those quantities enter into the 
dynamic equation and stress the fact that they should not be excluded a 



priori, particularly during the most rapid phases of the collapse. 

From (60) it is apparent that the production of density inhomogeneities 
is related to a quantity involving dissipative fluxes and shear viscosity. Thus 
if following Penrose we adopt the point of view that self-gravitating sys- 
tems evolve in the sense of increasing of density inhomogeneity, then any 
alternative definition for an arrow of time should include those variables. 

Finally, it is worth mentioning that we have considered the fluid to be 
neutral. The reason for this is that, as it can be easily verified, there is 
not terms coupling electromagnetic and dissipative variables, in the relevant 
equations. Therefore the role of electric charge in the dynamics of collapse 
is the same already discussed in [32]. 
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